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Abstract 

The encapsulation of hydrogen-terminated nanosized diamond fragments, the so-called diamon- 
doids into armchair single walled carbon nanotubes with diameters in the range of 1.0 up to 2.2 
nm has been investigated using classical molecular dynamics simulations. Diameter dependent 
molecular ordered phases were found for the encapsulation of adamantane (C10H16), diamantane 
(C14H20), and dihydroxy diamantane (C14H20O2). The same types of chiral ordered phases (dou- 
ble, triple, 4- and 5-stranded helices) observed for the encapsulation of Ceo molecules were also 
observed for diamondoids. On the other hand, some achiral phases comprising layered structures 
were not observed. Our results also indicate that the modification of diamantane through func- 
tionalization with hydroxyl groups can lead to an enhancement of the packing of molecules inside 
the nanotubes compared to non-functionalized compounds. Comparisons to hard-sphere models 
are also presented revealing differences, specially when more asymmetrical diamondoids are consid- 
ered. For larger structures (adamantane tetramers) we have not observed long-range ordering for 
nanotubes with diameters in the range of 1.49 to 2.17 nm but only a tendency to form incomplete 
helical structures. 

PACS numbers: 



I. INTRODUCTION 



The discovery of the Cqq fullerene^, carbon nanotubes^ and the first realization of encap- 
sulation of Cgo in carbon nanotubes^ have generated high interest in many research areas in 
nanotechnology. The large amount of experimental and theoretical works for these materials 
have been demonstrated their variety of functionalities making them a new and interesting 
class of nanostructured materials. 

The structure of single walled carbon nanotubes (SWCNTs) containing a linear arrange- 
ment of C 60 fullerenes, generically called "nano-peapods"^, has motivated experimental stud- 
ies of ordered phases of C 60 inside boron nitride^ and carbon nanotubes^ of large diameters. 
A model of cylindrically confined hard spheres has been employed to demonstrate the emer- 
gence of ordered phase formation due to packing maximization^. A systematic theoretical 
study using a continuum model was applied to the packing of Ceo inside SWCNTs and 
demonstrated that the formation of different ordered phases of Cqq molecules, resulting 
from the configuration of minimum energy of the system, is SWCNT diameter dependentP. 
These phases consist of zigzag, double and triple helices, and more complex arrangements. 
Some of these predicted ordered phases have been experimentally observed by Khlobystov 
et al. 6 during the encapsulation of Ceo fullerenes inside double walled carbon nanotubes 
with internal diameters of from 1.1 up to 2.6 nm. Atomistic studies using molecular dy- 
namics simulations of the encapsulation of C$$^ and smaller fullerenes (C20 and C2s)^into 
SWCNTs have shown a good agreement with the hard-sphere and continuum models. In 
addition, atomistic studies of the encapsulation of C70 and C78 also showed the dependence 
on the tube diameter, revealing similar types of ordered arrangements as the ones observed 
in C20, C28, and Ceo cases, despite of the asymmetrical shapes of C70 and C78 molecules^. 

Recently, other carbon nanostructures, known as diamondoids, have also attracted the 
attention of the scientific communitj^^. Diamondoids are essentially hydrogen-terminated 
nanosized diamond fragments (Fig. 1). Adamantane (Fig. 1(a)), diamantane (Fig. 1(b)), 
and triamantane, present in petroleum crude oil at low concentrations^, are the smaller 
diamondoids molecules, each one having only one isomer. The first structural elucidation 
of the smallest diamondoid, adamantane, occurred in 191Sp^. It was found and isolated 
from Czechoslovakian petroleum in WSS 1 ^^ but only in 1957 it was synthesized^. Small 
(1-2 nm in diameter) and large (> 2 nm in diameter) diamondoids are found in multi- 



2 



pies of adamantane. Larger diamondoids can appear in many shapes - rods, disks and 
even screws. Adamantane oligomers^ and polymers present high stability, strength, and 
hardness. Due to their flexible chemistry, pure diamondoids are particularly attractive for 
preparation of pharmacological active compounds with pronounced membrane activity due 
to lipophilicity of the aliphatic diamondoid core structur e 1 12 1 18 1 . On the other hand, function- 
alized diamondoids^ are ideal candidates for molecular building blocks in nanotechnologj^ni, 
and for molecular electronic applications when the incorporation of functionalized groups 
allows their attachment to metal surfaces to form self-assembled monolayers^. 

Mcintosh et al^ have studied the structural, electronic, and encapsulation properties 
of adamantane and diamantane using density functional theory. The encapsulation was 
investigated in (n,n) armchair SWCNTs ranging from n = 4 up to 8. They found that the 
encapsulation of diamantane into armchair SWCNTs with diameter larger than 9.5 A (n > 
7) may occur spontaneously without any energetic cost. As pointed out by Mcintosh et al^, 
carbon nanotubes are promising structures to be used in assembling diamondoids. However, 
investigations of the diamondoid encapsulation in larger diameter SWCNTs (n > 8) have 
not been yet carried out. For this reason, in this work, we have carried out atomistic 
molecular dynamics simulations to study the encapsulation of small diamondoids inside 
armchair SWCNTs to investigate the possibility of formation of ordered phases. 

A good agreement between the ordered phases predicted by the hard-sphere and con- 
tinuum models, and atomistic simulations has been observed for the C 60 molecule due its 
spherical symmetry^. On the other hand, one may not expect such agreement for asym- 
metric diamondoid structures. Aspects such as functionalization of diamondoids either by 
atom substitution or insertion, can lead to asymmetrical characteristics in geometrically 
symmetrical unmodified diamondoids. The asymmetry is more difficult to be addressed 
with analytical/continuum methods but can be easily and directly addressed in atomistic 
simulations. 

II. METHODOLOGY 

We performed classical molecular dynamics simulations using the universal force fiel d 23 * 24 * 
implemented in the Cerius 2 package^ which includes van der Waals, bond stretch, bond 
angle bend, inversion and torsional rotation terms. This methodology has been proved very 
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effective in the dynamical aspects of carbon based nanostructures ^ 27 * 28 l The dynamics of 
diamondoids molecules within a cylindrically confined environment produced by the SWCNT 
is mainly due to van der Waals interactions. A Lennard- Jones 6-12 potential was used to 
describe the van der Waals interaction between carbon, hydrogen and oxygen atoms, i.e., 
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where x is the separation distance between the atoms i and j, and 



(1) 



Aj = V D i D i x a = (V2)(^ + xj), (2) 

with i,j = C (carbon), H (hydrogen), or O (oxygen). The following parameter values were 
used™ D c = 0.105 kcal/mol, x c = 3.851 A, D H = 0.044 kcal/mol, x H = 2.886 A, D = 
0.060 kcal/mol, and xo = 3.5 A. 

We have considered (n,n) armchair SWCNTs with 8 < n < 16, and length of 12 to 15 
nm. Based on results for the CgcF, where the nanotube chirality presents a minor effect 
on the types of ordered phases, we also expect the same behavior for the encapsulation 
of adamantane and diamantane. The tubes were kept fixed during the simulations. We 
used the (8,8) nanotube as the smallest nanotube since the entry process of diamantane 
is energetically favorable for this nanotube 2 ^. We have investigated the encapsulation of 
adamantane T d (Ci H 16 , Fig. 1(a)), diamantane C 2 h (C14H20, Fig. 1(b)), functionalized 
dihydroxy diamantane (C14H20O2, Fig. 1 (c))^, and 1,3-adamantane tetramer (C 4 oH 58 , Fig. 

1 (d)jni. 

For SWCNTs with n < 12 we started the simulation by introducing the diamondoids 
one by one, until the complete filling of the nanotube. For the other cases, two or three 
linear arrangements of diamondoids have been initially inserted. The molecular dynamics 
simulations were then carried out in the canonical ensemble during hundreds of ps (time 
step=0.5 fs) at 300 K. For controlling temperature the Nose thermostat)^ was used. To 
avoid local minima, the temperature has been gradually increased until 700 K by steps of 
100 K per 50 ps. After that, a temperature decreasing of the system is carried out back 
to 300 K. At this final stage, all structures were geometrically optimized with convergence 
tolerances for energy and force of 10~ 4 kcal/mol and 5xl0~ 3 kcal/mol/A, respectively. The 
resulting optimized structures were used to identify the ordered phases. 
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(a) (b) (c) 




FIG. 1: Molecular structure of (a) adamantane, (b) diamantane, (c) dihydroxy diamantane, and (d) 
1,3-adamantane tetramer. Carbon atoms are represented in gray, hydrogen in white, and oxygen 
in red. 

III. RESULTS 

Adamantane and diamantane inside armchair SWCNTs with n = 8 and n = 9 were 
predicted to form a linear arrangements. Different ordered phases were obtained when n > 
9 to n = 16. The obtained phases for these diamondoids molecules are listed in Table 1. 
We can see that the same type of ordered phases observed for C 60 molecules can also be 
present in the encapsulation of adamantane and diamantane, i.e, zigzag, double, triple, 4- 
and 5-stranded helices. However, some Cgo achiral phases such as two- and three-molecule 
layers were not observed for adamantane and diamantane for the SWCNTs considered here. 
For the Ceo encapsulation^, these phases appear between double and triple helix phases, and 
triple and 4-stranded helix phases, respectively. For adamantane and diamantane, we can 
see a direct transition from double to triple helix phase and from triple to 4-stranded helix 
phase without these intermediate achiral phases. Excluding the linear and zigzag achiral 
phases, our simulation results suggest that adamantane and diamantane show a tendency 
to preferentially form chiral phases, at least until diameters up to about 2.2 nm. It is an 
indication that not only the nanotube diameter but also the structure of the encapsulated 
molecule play important roles on the ordered phase type that eventually is formed within 
the SWCNT. 

Figures 2 and 3 show some of the ordered phases obtained for adamantane and diaman- 
tane, respectively. For the (16,16) tube (Fig. 2), the adamantane molecules form a cylindri- 
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cal shell consisting of a 5-stranded helix on the internal nanotube wall. This phase has been 
predicted by Hodak and Girifalco as the minimum energy configuration of Cgo fullerenes 
inside (21,21) SWCNT 21 . It is important to notice the presence of several holes (Figs. 2 and 
3), corresponding to missing molecules, in the packing configurations. These holes along 
helix configurations appear to be most frequently observed for diamantane arrangements. 
Despite the presence of such defects, which cause a local disorder in the arrangement, it 
is still possible to observe some degree of organization that indicates the existence of an 
ordered phase. From the geometrical point of view, an analogy between adamantane and 
Cqo, and between diamantane and C70 can be made due to the more pronounced asymmetry 
presented by diamantane and C70 molecules. In fact, more disorder is observed in the ar- 
rangement presented by the encapsulation of diamantane than in the case of adamantante, 
similarly to the cases of ordered phases of Ceo and Cftf". This disorder plays an important 
role in the phase formation as seen for the specific case of (16,16) SWCNT encapsulated with 
diamantane (Figure 3). In this case, an incomplete 4-stranded helix is observed, acting as a 
intermediate phase between a 4- and 5-stranded helix, the former appearing in the (15,15) 
and the latter expected to be possibly formed in the (17,17) SWCNT. 

TABLE I: Packing phases obtained from molecular dynamics simulations of diamondoids inside 
armchair SWCNTs of different diameters. The arrangements presented by C60 molecules are also 
shown for comparison 9 . 



Nanotube Cgo fullerene 


Adamantane 


Diamantane 


Dihydroxy-diamantane 


(8,8) 




Linear 


Linear 


Linear 


(9,9) 




Linear 


Linear 


Linear 


(10,10) 


Linear 


Zigzag 


Zigzag 


Linear 


(11,11) 


Zigzag 


Zigzag 


Zigzag 


Zigzag 


(12,12) 


Zigzag 


Double helix 


Zigzag 


Double helix 


(13,13) 


Zigzag 


Triple helix 


Double helix 


Double helix 


(14,14) 


Zigzag 4- 


-stranded helix 


Triple helix 


Triple helix 


(15,15) 


Zigzag 5- 


-stranded helix 


4— stranded helix 


4— stranded helix 


(16,16) 


Double helix 5- 


-stranded helix Incomplete 4— stranded helix 


5— stranded helix 



In addition to pharmacological applications of functionalized diamondoids, Mcintosh et 
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(14,14) 




(16,16) 



FIG. 2: Double, triple, 4-stranded, and 5-stranded helix molecular phases obtained from molec- 
ular dynamics simulations of adamantane encapsulated in (12,12), (13,13), (14,14), and (16,16) 
SWCNTs, respectively. Different colors are used to represent the helical structures. 

al. have proposed that they can be used to increase the interaction between adjacent di- 
amondoids within a SWCNT^l The functionalization would enhance bonding and would 
allow the connection between encapsulated diamondoids, favoring a diamondoid polymer- 
ization inside the encapsulating tube. They proposed the substitution of terminal carbon 
atoms by boron and nitrogen in adjacent diamantanes in a linear arrangement. For this 
case, they obtained a value of 1.64 eV for the binding energy of functionalized diamantanes 
in contrast to 0.13 eV of the unmodified ones. While this type of functionalization has not 
yet been produced, other types have been already synthesized. An example is the dihydroxy 
diamantane (Fig. 1 (c))r™ Preliminaries calculations involving the reaction of two dihy- 
droxy diamantane molecules to form the CiztHigO dimer through the release of water have 
indicated a highly endothermic process, which would forbide the diamantane polymeriza- 
tion. In despite of that, we have investigated the encapsulation properties of the dihydroxy 
diamantane since hydroxyl group usually enhances the molecular reactivity that could allow 
diamondoid polymerization inside the nanotube through alternative reaction pathways. 

We have also observed the formation of ordered phases for the dihydroxy diamantane but 
with some differences with respect to the unmodified diamantane. The obtained phases are 
also shown in Table 1 and some of them presented in Figure 4. The dihydroxy diamantane 
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(14,14) 




(16,16) 



FIG. 3: Double, triple, 4-stranded and incomplete 4-stranded helix molecular phases obtained 
from molecular dynamics simulations of diamantane encapsulated in (13,13), (14,14), (15,15), and 
(16,16) SWCNTs, respectively. Circles indicate molecules that were found to not belong to any of 
the identified helices of the arrangement. 

forms a linear structure inside armchair nanotubes with 8 < n < 10. For n > 9 local 
disorder was observed to be more pronounced. A zigzag arrangement was observed when the 
molecules are confined into (11,11) SWCNT. The double helix arrangement was obtained 
for (12,12) SWCNTs. For the same SWCNT a zigzag pattern was found instead for the 
unmodified diamantane. We can also see a double helix arrangement in (13,13) SWCNT 
encapsulated with dihydroxy diamantane but with a different pitch of the one of the helices 
observed in (12,12) SWCNT. However, a similar pitch is found for the helices presented 
in (13,13) SWCNT encapsulated with dihydroxy diamantane and diamantane, respectively. 
Increasing even further the SWCNT diameter, complex helical arrangements are also formed 
such as a 5-stranded helix in (16,16) SWCNT, instead of the 4-stranded helix observed for 
unmodified diamantane molecules encapsulated into the same SWCNT. From these results 
we can conclude that the functionalization with hydroxyl groups in diamantane (Fig. 1 (c)) 
seems to optimize the packing. 

For the more asymmetrical diamondoid-based compound investigated here, i.e., 1,3- 
adamantane tetramer (Fig. 1 (d)) we have found that the smallest tube where it can be 
encapsulated is the (10,10) (Fig. 5). We have observed the formation of cylindrical shells 
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(12,12) 




(16,16) 



FIG. 4: Double, 4-stranded and 5-stranded helix molecular phases obtained from molecular dynam- 
ics simulations of dihydroxy diamantane encapsulated in (12,12) and (13,13), (15,15), and (16,16) 
SWCNTs, respectively. 

around of the SWCNT internal wall with some tendency to form helices for (11,11) up to 
(16,16) SWCNTs. However, in general, we could not see the presence of totally complete 
ordered phases during the simulation time used in our simulations. For instance, in Figure 
5, we can see an incomplete 4-stranded helix arrangement in (16,16) SWCNT where the 
incomplete helix (yellow) appears due to some misalignment of some molecules inside the 
nanotube (inside the circle). 

Analyzing the volume fraction^ as a function of nanotube diameter was possible to com- 
pare the Ceo packing with results obtained from the hard-sphere modeP. A good agreement 
with the hard-sphere model was obtained for SWCNTs with n > 15, when the ordered phases 
are three-dimensionaP. In order to make comparisons between the packing of diamondoids 
predicted by atomistic simulations and by the hard-sphere model, we have used the value 
of the available SWCNT diameter D = D n — A, where D n is the diameter of the SWCNT 
optimized with the universal force fielcF', and A = 3.38 A is the thickness of the nanotube 
wall. For the hard-sphere model D is written in units of the ball diameter^. For each di- 
amondoid curve we have expressed D in units of the its molecular diameter d assuming a 
spherical shape, i.e., d = 2(3\^/47r) 1 / 3 , where V is the molecular volume comprised by the 
Connolly surface 32 . The obtained volume values for adamantane, diamantane, dihydroxy 
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FIG. 5: Linear chain and incomplete 4-stranded helix molecular phases obtained from molecular 
dynamics simulations of 1,3-adamantane tetramer encapsulated in (10,10) and (16,16) SWCNTs, 
respectively. The circle indicates molecules that are misaligned with the helix pattern. 

diamantane, and Cgo were (in A 3 ) 156.8, 205.1, 217.1, and 490.6, respectively. With the 
volume calculated for the C 60 molecule by this method, the corresponding value of the C 60 
diameter is in agreement with the intermolecular distance 0.97 ± 0.02 nm obtained for the 
linear arrang ement in the (10,10) SWCNT®. 

Figure 6 shows the volume fraction as function of D for the encapsulation of adamantane, 
diamantane, dihydroxy diamantane, and C§q molecules. We can see that all diamondoids 
for the linear arrangement in the (8,8) SWCNT (D = 1) present a volume fraction smaller 
than 2/3 (the expected value for the encapsulation of hard spheres in a hard cylinder of 
the same sphere diameter). Due to the spherical symmetry presented by the Ceo molecule, 
a reasonable agreement with the hard-sphere model is seen for D > 1.7 . For D < 1.7 an 
increase in the volume fraction is observed for the Cgo packing. In this region the phases are 
basically two-dimensional and, therefore, the number of encapsulated fullerenes is limited 
by the nanotube length. As a direct consequence, increasing the nanotube diameter causes 
the volume fraction value to decrease. The behavior of the volume fraction shown by the 
adamantane encapsulation is qualitatively similar to the Ceo, which can be associated with 
an approximated spherical shape also shown by the adamantane. In both cases we can 
see an increase of the volume fraction during the transition of the zigzag to double helix 
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FIG. 6: Volume fraction as function of the available diameter D expressed in units of each molecular 
diameter. The lines connecting the points are only to help visualization. 



arrangements (D ~ 1.73). Furthermore, the behavior for the adamantane is qualitatively 
different for the other diamondoids considered here. For both diamantane and dihydroxy 
diamantane, the volume fraction begins to increase after D ~ 1.93 where the double helix 
arrangement is present ((13,13) SWCNT). The effect of functionalization also reflects on 
the volume fraction values. For the functionalized diamantane (dihydroxy diamantane) we 
can see higher values of the volume fraction than the unmodified diamantane, indicating a 
packing enhancement produced by the introduction of hydroxyl groups in diamantane. This 
enhancement is more evident for the achiral linear and zigzag phases (D < 1.7). 



IV. CONCLUSIONS 



Using atomistic molecular dynamics simulations we were able to identify molecular or- 
dered phases on the encapsulation of diamondoids into armchair single walled carbon nan- 
otubes. We have verified that the arrangement types depend on the diameter of the encap- 
sulating nanotube as well as the structure of the encapsulated molecule. For the SWCNTs 
considered in this work the investigated diamondoids have not shown some achiral phases 
that have been previously predicted on the encapsulation of Cgo molecules. The asymmetry 
presented by the diamondoids can lead to local disorders on the arrangement that do not 
affect the long-range ordering. An enhancement on the packing and a favoring to chiral 
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phases were obtained in the encapsulation of functionalized diamondoids. The more packed 
phases shown by functionalized diamondoids may serve as candidates for diamondoid poly- 
merization inside carbon nanotubes. The encapsulation of larger structures (1,3-adamantane 
tetramers) was only found for tubes with diameters larger than 1.36 nm ((10,10) SWCNT). 
For those diamondoid-based polymers we have not observed long-range ordering for diam- 
eters larger than 1.49 nm but only a tendency to form incomplete helical structures. This 
long-range ordering is precluded by misaligned molecules within the helical arrangements. 

We acknowledge the financial support from the Instituto de Nanosciencias, THEO- 
NANO, IMMP/CNPq, the Brazilian agencies FAPESP, Capes, and CNPq. 
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